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Q . Abstract 

We consider two loop sunset diagrams with two mass scales m and M at the threshold and 
pseudotreshold that cannot be treated by earlier published formula. The complete reduction 
to master integrals is given. The master integrals are evaluated as series in ratio m/M and in 
e with the help of differential equation method. The rules of asymptotic expansion in the case 
when q 2 is at the (pseudo)threshold are given. 

r- 1 Introduction 

o 

The sunset diagram plays a key role in the 2-loop calculations with masses. Despite the fact 
that a lot of investigation has been devoted to the sunset diagram there still remain drawbacks. 
In Ref. Q general reduction procedure is given in the case, when external momentum q and 
internal masses are arbitrary. But in the case when q 2 is equal to the threshold or one of the 
pseudothresholds values the formula of [Q] become unapplicable. Therefore we turn to paper 
i-^h where the reduction was given specifically for the (pseudo)threshold kinematics. However 

in two cases shown in Fig.l the reduction of Q fails. 

In present paper we consider calculation of these two special cases. These integrals naturally 
arise in the threshold problems with given mass hierarchy. An immediate typical example, 



C3 



where it is the case, is a problem of matching vector and axial QCD currents to NRQCD ones 
with two heavy quark mass scales m <C M. Another example is the calculation of masses of 
the heavy gauge bosons in 2-loop approximation §, when q 2 is equal to m 2 z or m 2 ^. 

For the calculation of master integrals there are mainly three methods: 1) direct evaluation 
using a or Feynman parameter representation, 2) solving master differential equation in exter- 
nal Mandelstam variables, which can be written for the master integrals of any Feynman graph, 
3) applying various asymptotic expansions @]. Here, we will demonstrate the strong and 
weak features of the each mentioned method on the example of our particular problem and 
will advocate that a certain mixing of these methods can give us a desired answer in the most 
easiest way. 

Now let us introduce notation, which will be used later in this paper and define master inte- 
grals we are going to calculate. For the two-loop sunset with arbitrary masses and propagator 
indices we have 

2 1 ff d d kd d l 
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where d = 4 — 2e is the dimension of space-time. 

We will discuss only two special cases of the sunset integrals which are shown in Fig.l: 

(a) mi = M, ni2 = 0, 1113 = m, q 2 = (M + m) 2 (threshold case), 

(b) mi = rn, ^2 = m, 1713 = M, q 2 = M 2 (pseudothreshold case). 

M 

>— q 

m ' M 

Case (a): cf = (m+M) 2 Case (b): cf = iVf 

Figure 1: Sunset diagrams 

With the help of recurrence relations []|], which will be given in next section, any integrals 
of the first type with arbitrary propagator indices can be reduced to two master integrals: 
Jm((m + M) 2 ) and Jn2((m + M) 2 ), whereas in the second case to two master integrals: 
Jin (M 2 ) and J 211 (M 2 ). 

Another problem is the evaluation of the master integrals themselves. The general repre- 
sentation for the sunset diagram was obtained in Ref. || in terms of hypergeometric Lauricella 
function. For the practical purposes however one needs the e-expansion of these formula which 
is not a trivial task. 

The result for the threshold and pseudothreshold values of the sunset integrals with three 
arbitrary masses has been ontained in Ref. |l0| up to 0(1) in e-expansion. From these results 
we can easily obtain the finite parts of our integrals. However we established that in the 
reduction procedure for the matching of QCD to NRQCD currents one has to know also 
O(e) part for integrals and 0{e 2 ) part for jw integrals. It is not easy to push forward 



the approach of [10] in order to evaluate these needed parts. Instead we use the differential 
equation method [pd|]-|13|| (see also |14[| ). We solve differential equations as series in m/M 
and the desired number of coefficients can be always obtained. In order to find the boundary 
conditions for the solutions one can use two methods: the representation of Ref. ||] can be 
expanded in the limit m/M — » and the asymptotic expansion procedure can be applied at 
the threshold and pseudothreshold. 



2 Recurrence relations at the threshold and pseu- 
dothreshold 

Here we give the recurrence relations § for the sunsets of two types introduced earlier. While 
the derivation of these recurrence relations is straightforward, they were not considered in 
the literature in detail till now. On the other hand they represent missing pieces to com- 
plete the generalized recurrence relations of Tarasov Q and threshold relations of Davydychev 
and Smirnov Q for the sunset diagrams in the case of threshold with one zero mass and 
pseudothreshold with two equal masses. Most of the formula below can be derived just by 
combining and reexpressing appropriate recurrence relations of |] and @. 
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2.1 JMOm 

In this case we can not directly apply recurrence relations for general masses obtained in [||. 
The point is as one of mass becomes zero the relations of |2| are degenerate. Instead we derive: 

{d - 2v 2 - 2) u 2 2 + = -2m 2 v 3 {v 3 + 1) 3++ + (d - 2v 3 - 2) u 3 3 + , (2) 

2M 2 v 1 {v 1 + 1) 1 ++ = 2m 2 v 3 {v 3 + 1) 3 ++ - {d-2v 3 - 2) v 3 3 + + (d-2v x - 2)^1+ , (3) 



2m 2 (M + m){3d - 2v x - 2v 3 - 7)u 3 {u 3 + 1)3 ++ = 

m((d -v 3 - 3){2d -v x - 2v 3 - 4) + {d - v x - u 3 - 2) {3d - 2v 2 -2v 3 -Q 

+M{2d 2v 3 - 5) (2d - v x - 2v 3 - 4)1 v 3 3 + 

+{d -v 3 - 2){m{d -v 3 -3)+ M{2d -v x - 2v 3 - 5)) v x \+ 

+2m 2 Mv 1 v 3 {v 3 + l)l+2"3 ++ 

+ \-m{d -v 3 -3)- M{2d -v x - 2v 3 - 5)1 v 1 v 3 l + 2-3 + , 



(4) 



2mM{m + M){3d - 11) J 2 i 2 = 

-{d - 3) (m{d - 3) + M{2d - 7)) J 112 - {d - 3) (m{2d - 7) + M{d - 3)) J 2U 
m + M 



Am 2 M 2 



{d-2) z {2d-7)J l01 . 



(5) 



Here as usual j + (or j ) means the operator raising (or diminishing) index on the j-th line. 

With the help of ©-(1) we reduce all integrals to three sunset diagrams Jm, J\\ 2 and J 2 \\ 
plus products of 1-loop tadpoles. In addition there is a relation between these three sunset 
integrals and one of them can be eliminated 



M{m + 2M) J 211 + m(M + 2m) J 112 
This finishes the reduction procedure. 



3d 



I {d - 2)2 7 

JlU + 4mM(d - 3) Jl01 



(6) 



2.2 JinniM 

Here we are faced with pseudothreshold problem and again in this case general recurrence 
relations by Tarasov [|J] and those for the threshold problems by Smirnov and Davydichev [|| 
become degenerate. The reduction procedure for this type of integrals has been studied in [16], 



where the topology with three and four lines has been considered. These integrals were also 
considered in ||, where asymptotic expansion has been used. For the completeness we give 
alternative reduction formula below. 
First, using relations 

2m 2 v x {v x + 1) 1 ++ = {d-2v x - 2)v 1 l + - {d-2u 3 - 2) v 3 3 + + 2M 2 v 3 {v 3 + 1)3 ++ , (7) 
2m 2 v 2 {v 2 + 1) 2 ++ = {d- 2v 2 - 2) v 2 2 + - {d - 2v 3 - 2) v 3 3 + + 2M 2 v 3 {v 3 + 1) 3++ , (8) 

we reduce indices of lines 1 and 2 to one or two. Thus there remain only integrals J\\ Vz , Ji2^ 3 
and J 22v3 . 
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For J\2u 3 we have 



4(AT - m z ) i/ 3 3 + = 
For Jrivz we have 

4m 2 (M 2 - m 2 )(d - v 3 - 3) = 
m 2 (d-v z - 3)3" + 



(d - 3) 2" J/ 3 3 + - 1+3" + d - 3i/ 3 



-3c? 2 + d(13 + 7i/ 3 ) 



-12 - 17z/ 3 - 4z/f J + M 2 (2d - i/ 3 - 6)(d - i/ 3 - 2) 
+ [-m 2 + M 2 (d-u 3 - 2)1 (d - 3) i/ 3 3+l~2~ 
+ [m 2 (d- 3) - M 2 {2d-v 3 - 6)1 i^3 + l~ . 



For Jii^g we have 



++ 



4M 2 (M 2 - m 2 )(d - z/ 3 - 3)i/ 3 (i/ 3 + 1) 3 

- m 2 (d - 2^3 - 3) (2d - 2u z - 5) 

+Af 2 (3d 2 - d(17 + 7i/ 3 ) + 24 + 19z^ 3 + 4z/f)l z^ 3 3+ 

+ [ - m 2 (d 2 - d(3 + 5^ 3 ) + ^(4^3 + 11)) -M 2 {d-u 3 - 2)u 3 

+ [- m 2 (d - 3) + Af 2 ^ 3 j z^3 + 2 + l~ + m 2 (d - u 3 - 3) 3" 1+2+ . 

Finaly there is relation between three integrals 



(9) 



(10) 



„ ,9 T 9 T 3d — 8 _ 

M 2 J 112 + m 2 J211 = — - — J 



111 



(d - 2) 2 
8m 2 (d-3) 



110 • 
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(12) 



This finishes the reduction. 



3 Master differential equation 

It is known, that a general two-loop sunset topology has four master integrals jjjj: one with 
unit powers of propagators and three with a dot placed on one of the lines. For master integrals 
considered in this paper this number is however smaller, which is due to a symmetry of the 
mass distribution on the lines. In addition in case of JMOm the index of massless line can 
be always reduced to 1 (see fij). It can be shown, that these master integrals satisfy a sytem 
of linear nonhomogeneous differential equations in q 2 (q being external momentum) |l3| or in 
masses [11]. However, in our case we want to write differential equations in r = m/M, where m 
and M are internal masses and q 2 lays at threshold or pseudothreshold. Differential equations 
of this type were considered earlier in fl 



3.1 JMOm 

In this case we have two master integrals (dm and dm) and the system of two 1st order 
equations. It is convinient to rescale integrals introducing dm and dm by 



d 



111 



M 



2d- 



T 2 (3 - d/2) dm and d m = M M " 8 r 2 (3 - d/2) d U2 . 



(13) 
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Then the differential equations read 

(r + i)( r + 2) — Jin - 6r(r + 1) J 112 -(d + 2{d - 3)r - 4J J m - _ _ - = , 

(r + l)(r + 2)^-J 112 + -(2r 2 - 5dr + 20r - 4d+ 14)J 112 ' ' " 

or r 

(d-3)(3d-8) ? 8r d - 5 
2r (<i-4) 2 

Using this system one can write the second order differential equation for Jm: 
(r + i)2 (r + 2 ) 2 |l J m _ H (r + 1)(r + 2) 2^ + {d _ 4)r _ 3 ) 9_j ni 

+-(d ~ 3) 2 (r + 2) 2 (d - 2r- 4) J m - — ^ (r + 2) V" 4 = 0. (15) 
r (a — 4j z 

As usual we search the solution of (|l5|) as a linear combination of two solutions of the homo- 
geneous equation and the solution of the nonhomogeneous equation. We will find the solution 
as a series in r, namely, in the following form 

^m=£r^f>«(d)A (16) 

i \n=0 ) 



By substituting (|16|) in Eq. fllq), we find for leading exponents three allowed values a\ = 
0, «2 = 2d — 5, corresponding to the two independent solutions of the associated homogeneous 
equation and «3 = d — 2, as required by the nonhomogeneous part of Eq. (H). So the 
only thing we need yet to do is to find the coeffients in front of two independent solutions of 
homogeneous part of Eq. (||) using boundary conditions. One equation for these coefficients 
can be obtained from the value of Jm master integral at r = 0, which can be written in 
terms of T-functions. One can not use Jm at r = as the second boundary condition. The 
reason is that the latter integral becomes infrared divergent and this divergency is regularized 
by e {d = A — 2e) and not by r as in the limit r — > 0. In order to obtain a second boundary 
condition we need explicit expansion of function Jm in r up to the third order, which can 
be obtained by analysing representation of this integral in terms of Appel function F4 or 
performing asymptotic expansion in small mass ratio r. Expansion of the Appel function will 
be considered in Appendix A, while asymptotic expansion technique is discussed in Section 
4.1. 

Having obtained expansion in r up to order 0(r 3 ) using one of the mentioned methods, 
one can easily reconstruct all other expansion coefficients from differential equation. Returning 
back to functions Jm and Jm instead of Jm and Jm we give the first seven coefficients of 
expansion: 

s^M-^J ni = ("-B)r^2r-5 _ 1 (n + 20<2 , + | 

-l(16L 2 - 48L - 12C 2 + 7)r 2 - ^(9L 2 + 15L + 5A( 2 - 8)r 3 + ^(24L 2 + 20L + 144( 2 - 3)r 4 

o y o 

— ^-(450L 2 + 255L + 2700C 2 - 16)r 5 + ^(360L 2 + 156 J + 2160C 2 - 5)r 6 
225 72 

+ j55_^_l|3 + l (20C2 + 11)r + AL3_ 6i2+2(C2 + 7)i 
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and 



+^(276C2 + 20% + 321)) r 2 + (4I 3 + H|! + (406 - ?f)L - ^ - 24( :1 + ^ j r 3 

+ (-61*-f + (^ -606)1 + 176 + 366- ff)r< 

+ L + + (806 - ^1)L - Mi _ 48 <3 + ^) r = 
^ 15 v s 450 ; 15 s 27000 y 

+ (-101* - ^ + ( ™ _ 100C2)i _ S + 60(3 _ ^ j r «) + 0(r , + 0(E , (17) 

d 

1 2L — - 1 

e 2€ 7BM 4ej ii2 = -— + 2 + ( g C2 + 2 ) + 4L _ 2L 2 + ( _ 2L 2 _ ^ _ ^ + 2 y 

+(3L 2 + 3L + 18C 2 - i)r 2 + (-4L 2 - ^ - 24 C2 + ^)r 3 + (5L 2 + ^ + 30C 2 - ^)r 4 
+( _ 6L 2 _ lp _ 36C2 + ^ )r 5 + (?L 2 + JL + _L (756C2 _ 1))r 6 

(1 4L 3 
— (66C2 + 52C 3 + 66) + 8L + 2LC2 - 4L 2 + — + (4L 3 + 4L 2 + 40C 2 L - 24L 

-8C2 - 24C 3 + 14)r + (-6L 3 - 5L 2 - 6OC2L + 37L - 6C2 + 36C 3 - —)r 2 

( - 22L 2 , r 440L 68C2 , 1 1277 \ , 
+ (^ + — + 8061 - — + -f - 486 + — ) r» 

+ (-101= - ^ _ 100(2 1 + i»™ - 42 C2 + 60 (s - 

+ f 12i 3 + + 120f2i _ 3186 _ J5179\ 

\ 5 s 75 5 w 1000 / 

/ Jr , 529L 2 i , r 6418L 1307C 2 _ 4823273 \ fi \ ^, 7n 2n 

+ -14L 3 140C 2 L + — + 84C 3 r 6 + O r 7 + O e 2 , 

^ 30 s 75 15 s 54000 J J v ; v h 

(18) 

where L = logr. The 0(1) part of Jm is in agreement with [1C]. 

3.2 JmmM 

This diagram was recently studied in Ref. (T^j by means of differential equation method in 
the regime when m 3> M, while we are interested in the case m <C M. Here we also have 
two master integrals and the system of two differential equations. Again introducing rescaled 
functions according to ( |l3|) we have 

d ~ 

— Jm - 4rJ 2 n = 0, 
or 



r 4 



r 5 



(19) 

2r z ) _ n 

dr" z±1 ' r " z±1 ' 4r (d — 4) 2 



, 2 ,J= r 2 (13 -4d) + 2d- 7 - (d-3)(3d-8) ? 2 r d - 7 (r d - 2r 2 ) 

( r ~~ 1) 7^-^211 H J211 H "j Jill 73 TTo = 0. 
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Using (|19|) the corresponding second order differential equation for Jm looks like 

, 2 . d 2 2(d-3)(2r 2 - 1) d . x 8r d-6( r d _ 2 2j 

( r " !)^2 J m " 7 J m + ( d " 3 )( 3d " 8) (d _ 4)2 ^ = 0. (20) 



To find a solution of this equation we again use Ansatz (|16|) for the most general form of 
solution at r — » 0. For leading exponents aj there are four allowed values a\ = 0, a 2 = 2<i — 5, 
corresponding to the two independent solutions of the associated homogeneous equation and 
«3 = d — 2, Q!4 = 2d — 4, as required by nonhomogeneous part of Eq. (p0|). All other steps in 
this case are in one to one correspondence with those considered in previous subsection. The 
first boundary condition is given by the value of master integral Jm at r = 0. In order to find 
the second boundary condition we need explicit expression for r-expansion of Jm up to the 
third order, which can be obtained from threshold large mass assymptotic expansion of Jm. 
As a result of all these steps we have the following expressions for our master integrals 

e^M- 2+4 7 m = + 4(4L ~ 3)r2 ~ 5 - i(H + 20C 2 ) - (4L 2 - 12L - 3C 2 + 5)r 2 

2e z 4e 8 

-\(8L 2 - 12L + 8C 2 + 7)r 4 - -(12L - ll)r 6 + e fbh UCs 25(2 



V ' 18 v ' \16 3 4 

+ - 12L 2 + 4(C 2 + 7)L + 9(2 + ^ - 3^J r 2 - 32C 2 r 3 

+ (4L 3 - 5L 2 - 7 ± + 4 C2 - 4C 3 + f ) r 4 + ^ + ^(8L - 9( 2 - 14)r 6 ) 

2 /949 55C 2 55C 3 1 a, a ■> 
+£ (-32" " -T " "IT - 720^ ^ 296r " ^ + 3 ° 3) 

+ {-^- + 8L 3 - 4(C 2 + 7)L 2 + ^(9( 2 - 2C 3 + 45)L + 31( 2 + 26C 3 + 19^ r 2 

64, , , / 14L 4 , 3 N 2 , 37, 
+_( 6L + 12 bg 2 - ll)C 2 r 3 + f — + 6L 3 + (- - 2C 2 )L 2 + (3C 2 + —)L 

+ 8C3 - ) r 4 - |(60L + 1201og2 - 77)C 2 r 5 
+— (288L 3 - 792L 2 - 108(2C 2 + 11) L + 1926C 2 - 1296C 3 + 5417)r 6 ) + 0(r 7 ) + 0(e 3 ) (21) 



and 



e 2 ^M 4 V 211 = -± + ^-1 + \ + ^ + AL - 2L 2 - (2L 2 -2L + 2( 2 + l)r 2 

+(4L 3 - 2L 2 - 6L + 4C 2 - 4( 3 + ll)r 2 + 8C 2 r 3 + \{2AL - 27( 2 - 38)r 4 + 

9 5 

. L t 83 . fi \ 2 /49 37C 2 37C 3 289vr 4 4 . 

or3 2/ 4 

-2L 2 C 2 - 8L 2 + — — + 48C 2 (2L + 4 log 2 - 3)r 



7 



+ 



' ' 1 ' ' L * - 2(C 2 - 3)L 2 + 2(C 2 + 5)L - 17( 2 + 8( 3 -^- 5s) r 2 



+ 



90 



■f (12L + 24>o g 2 - 13) f2 ^ + - 3L= - ( 6 + f )I + ^ _ 66 + 



- — (120L + 2401og2- 199)C 2 r 5 

/4L 3 5L 2 1 , . 37C 2 237511 \ 

+ 150G + 203 L + — - - 2C 3 + 

\ 9 9 450 v s ; 180 s 81000 I 



r 6 + 0(r 7 ) + 0(e 3 ) , (22) 



where L = logr. The 0(1) part of Jm is in agreement with [1C]. 



4 Asymptotic large mass expansion at the threshold 

In this section we consider asymtotic large mass expansion for master integrals introduced 
above. A type of expansion one needs to perform in order to obtain an analytic expression for 
the master integral of case (b) was already considered in Ref. [18] and |J and one may just 
follow along the lines of procedure described there. However, in the case of master integral 
JMOm a somewhat different prescription for setting loop momenta is requiredP] and thus we 
will consider this case in detail below. 

In order to establish the expansion procedure we use " the strategy of regions" p^] . Let us 



remind you first a general prescription for large mass expansion at the threshold of Ref. [18]. 
We consider a general case of threshold Feynman integral Fp, corresponding to a graph F when 
the masses Mj and external momenta Qi are considered large with respect to small masses 
rrii and external momenta qi. We are interested in a case, when external momenta Qi are on 
the following mass shell: Q 2 = (J2j a ijMj + J2k bik m k) 2 - Here dy and bik are some numbers. 



It is just the generalization of the on-mass-shell condition of Ref. [18]. Then the asymptotic 
expansion in the limit Qi, Mi — ► oo takes the following explicit form M 



F r (Qi,Mi, q umi \ e) Ms ~°° £ M^F r (Qi, Afi, ft, m*; e). (23) 



,W, — x 

7 

Here the sum runs over subgraphs 7 of T such that 

(a) in 7 there is a path between any pair of external vertices associated with the large external 
momenta Qf, 

(b) 7 containes all the lines with the large masses; 

(c) every connectivity component 7,- of the graph 7 obtained from 7 by collapsing all the 
external vertices with the large external momenta to apoint is 1PI with respect to the 
lines with small masses. 



Operator 7W 7 in ( [23]) can be written as a product UiAi^ over different connectivity com- 
ponents, where are operators of Taylor expansion in certain momenta and masses. In 
what follows we will distinguish the connectivity component 70, which is defined to contain 
all external vertices with large momenta. For connectivity components 7«, different from 70, 



2 Asymptotic expansions in momentum space are not invariant under the redefinition of loop momenta contrary 
to the expansions, performed in a-representation for Feynman integrals. Therefore special care should be taken 
choosing a correct one set of momenta. 
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the corresponding operator performs Taylor expansion of the Feynman integral in its 
small masses and external momenta. To describe the action of A4 70 one uses representation of 
7o-component in terms of a union of its 1PI components and cut heavy lines (that is subgraph 
becomes disconnected after a cut line is removed). Here we can again factorize «M-y and the 
Taylor expansion of the 1PI components of 70 is performed as in a case of other connectivity 
components 7,. As for the action of operator M on the cut lines there are two different cases. 
Let P + k be the momentum of a line with large mass Mj, where P is a linear combination 
of large external momenta and A; is a linear combination of loop momenta and small external 
momenta. Then the mentioned two cases can be written as follows: 

• P 2 = Mf and A4 for this line is given by 



M = % ' 



(24) 

x=l 



zk 2 + 2Pk 

Here T x denotes the operator of Taylor expansion in x around x = 0. 

• P 2 / M 2 and the operator Ai reduces to the ordinary Taylor expansion in small (with 
respect to this line) external momenta. 

As for the optimal set of internal loop momenta we choose a rule when a large external momenta 
are divided between lines with masses Mj and in order to satisfy the following conditions: 
P 2 = M 2 and P? = mi. We do not know whether such separation is always possible or not, 
but in most cases of interest it certainly works. As an example in the next subsection we will 
consider an expansion for master integral Jm((m + M) 2 ) from our case (a). 

4.1 Second boundary condition for JMOm 

For master integral Jm, according to our prescription for choosing internal momenta we have 
the following expression: 

1 rr d d kd d l 
111 [k 2 }[(k + l) 2 + 2a(k + l,q)}[l 2 + 2blqy (25) 

Here q is an external momentum; a = M/(M + m) and b = —m/(M + m). From the general 
formula ( |23|) in the asymptotic expansion of master integral Jm in the limit m/M — > we 
have four subgraphs^]: 

1) graph T itself; 

2) subgraph 71 consisting from lines with masses M and m; 

3) subgraph 72 consisting from lines with mases M and zero; 

4) subgraph 73 consisting from one heavy line. 

For 1) we expand the integrand around m = with the result 

v MH^r rr d d kd d i{i • or , , 

^ TT d J J [k 2 ][(k + l) 2 + 2a(k + l,q)][l 2 } n + 1 - 



3 Here T is the graph, corresponding to Feynman integral Jm. 
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Each term in this expansion can be evaluated rewriting it via scalar integrals with shifted 
space-time dimension jD] and for first three of them we have: 



1 ff d d kd d l(lq) 



■2 



a q ~JVi+i, 



TT d J J [k 2 ] [l 2 ] 2 [(k + l) 2 + 2a(k + l,q)] 

l_ rr d d kd d l(lg) 2 _ d+2 _ £ d+2 2 4 d+4 

n d JJ [k 2 ][l 2 ]3[(k + l) 2 + 2a(k + l,q)]- 2^ 231 2 J V 132 + 4a J V 333 , 

J_ ff d d kd d l{lqf _ 3 6 TV d+6 , Q „4 Tv d+4 , o„„4-r V d+4 

n d JJ [k 2 ][Pn(k + l) 2 + 2a(k + l,q)] ~ Q 444 + Q 243 + Q 342 ' 

where 

JV ^ 3 w dJJ [k 2 }vi[l 2 ]^l(k + l) 2 + 2(k + l,q)}»3 q 

The latter can be easily expressed in terms of T- functions for arbitrary values of u\, v 2 and v^. 
The Taylor expansion for subgraph 72 gives 

^(-1)" ff d d kd d l(k 2 + 2M + 2akqr = , , 

^ Q 7r d J J [k 2 ][l 2 + 2alq] n + 1 [l 2 + 2blq] [ ' 

In case of subgraph 73 we have the following expression 

^ (-1)" ff d d kd d l(l 2 +2kl + 2alq) n 

^ TT d J J [k 2 ] [k 2 + 2akq] n+1 [I 2 + 2blq] ' [ ' 

Each term of this sum is just a product of two one-loop integrals and hence can be easily 
evaluated. 

An expansion for subgraph 74 leads to the following result 

v tDl ff d d kd d i((k + i) 2 r , , 

^ 7T d J J [k 2 }[l 2 + 2blq}[2a(k + l,q)] n + 1 ' 



Here we can see, that it is the most difficult from computation point of view type of con- 
tribution, as we are dealing in this case with eikonal integrals. As one can easily see, the 
evaluation of each term from this contribution can be reduced to the evaluation of integrals of 
the following type 

1 ff d d kd d l(kl) m 

(30) 



ir d JJ [k 2 ][l 2 +Zblq][(k + l,q)] n ' 

In order to calculate these integrals it is naturally to express the products (kl) m in terms 
of traceless products (kl)( m ^ = k( a ' m Hr am -ffl. Then we notice, that integration over k of 
traceless products fc( a > m ) times the part of integrand, which depends only on k results in 
expression propotional to the traceless product q( a ' m \ Thus, we can replace the factor (kl)^ 
by (qk)( m ^ (ql)( m ^ / (qq)( m ^ and finally replace the factors involved through ordinary products qk 
and Iq. After performing these steps we reduce the evaluation of integrals (j30|) to the integrals: 

1 ff d d kd d l(lq) m 

(31) 



ir d JJ [k 2 ][l 2 + 2blq][(k + l,q)] n ' 



Here a is a collective index representing m Lorenz indices ai, oez, . . . , a r , 
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These integrals with the use of integration by parts identities (d— 3 — n)n + + {n+l)rn + n ++ = 0, 
where n(m) + are operators, which increase corresponding indices n or m) can be further 
reduced to integrals of the form 

J_ rr d d kd d l(lq) m , , 

ir d JJ [k 2 ][l 2 + 2blq][(k + l,q)Y [ ' 

To take these last integrals one can first perform ko and Iq integrations using Cauchy theorem 
and remained angular integrations are trivial since there are no products kl in the integrand 
left. 

However, as our calculations showed, it is far more simple just to calculate using asymtotic 
expansion this master integral up to the order (m/M) 3 and then use this result as the nessesary 
boundary conditions for the solution of master differential equation. It takes then much less 
efforts to construct the expansion of this master integral up to any order in r = m/M using 
differential equation method. 



5 Conclusion 

We considered two loop sunset diagrams with two mass scales m and M at the threshold 
and pseudotreshold that cannot be threated by earlier published formula |l|, ||]. The complete 
reduction to the master integrals is given. The master integrals are evaluated as series in ratio 
m/M and up to needed in applications order in e with the help of differential equation method. 
The rules of asymptotic expansion in the case when q 2 is at the (pseudo)threshold are given. 

Acknowledgements. We would like to thank K. Chetyrkin, A. Davydychev, VA. Smirnov, 
M.Yu. Kalmykov and O.V. Tarasov for fruitfull discussions of the topics discussed in this paper 
and valuable comments. Authors acknowlege Universitat Karlsruhe for the warm hospitality. 



A Expansion of F4 Appel function: second boundary 
condition for JMOm 

In order to obtain the second boundary integral in case of JMOm we can use the representation 
of the sunset diagram in term of Lauricella function which can be obtained from || . However 
since one mass is zero, the Lauricella function simplifies to Appel function: 

-r(i- e) r(-i + ,>r(-i + *)* + ; (33) 

and in our case q 2 = (m + M) 2 . We need the expansion of this integal up to the order 0(r 3 ) 
and 0(e). Then the rest coefficients can be found from the differential equation. 



In order to expand (33) in series over m/M we use the following representation for F4 

i)k{b)k x k 
(c) k kl 



, a, b \ ^ (a) k (b) k x k fa + k,b + k • 



k=0 
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In (^) function 2 -Pi has to be transformed to argument 1 — y. Then we have rot the r.h.s. 

(r = m/M) 



„2-2£-p2 



r 2 (-i + £ )£ 

k=0 



(l)k(e) k r 2k 
(2-e) k k\ ' 



r (x^r^-O 2^1 ( J tV 2 f ; -K2 + r)) 
+ [-r(2 + r)] 1 - 22 -^ r ( 2 r; fc ™) 2^1 (^r-t-tr' ; -K2 + r) 

-r(-i + £ )r(i - e )r(-i + &) £(-i + 2 £ ) fe — 

k=o K - I 

F I 2-e,3-4e-2fc \ Ei / -l+2e+k, e+k . _ r (0 _|_ r \\ 
1 \^3-3s-fc,2-2e-fey 2 2e+2fe ' ' \^ > ' ) I 

+ [-r(2 + r)] 3 " 4£ - 2fc r (^U+t^) 2^1 (^I^r* ; -K2 + r)) } . (35) 



The first and third terms in fl35| ) can be easily expanded since the series in k and 2 -Pi series 
can be truncated at the given order. In the second and fourth terms we still can truncate 2 -Pi 
series, however, we cannot truncate the sum over k because of factor [— r(2 + r)]~ 2k . Thus we 
have to resum the whole fc-sum. Thus in order 0(r 3 ) we have for the second term 



-r)- 2£ (+r)- 2£ 



2v^F 

and the fourth term 

(-r)~ 4e 



r(2 - e)r 2 (-l + e)r(-l/2 + e) 2 F l ~ k ; 1 



20Fr(g) r(2 " £)r2(1 " e)r(_1 + e)r( " 3/2 + 2e)r(_1 + 2e) 2i?1 ( 



3-3e-fc, 2-2e — k . 



4-4e-2fc 



;i ■ 



Each of this two terms has an imaginary part but it cancels in the sum of the two. Adding 
contributions from the first and third terms of (j35|) and expanding in s we get 0(r 3 ) term of 
the formula (p7|). 
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